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Abstract 

The laws of quantum-critical scaling theory of quantum hdelity, dependent on the 
underlying system dimensionality D, have so far been verified in exactly solvable ID 
models, belonging to or equivalent to interacting, quadratic (quasifree), spinless or 
spinfull, lattice-fermion models. The obtained results are so appealing that in quest 
for correlation lengths and associated universal critical indices ly, which characterize the 
divergence of correlation lengths on approaching critical points, one might be inclined 
to substitute the hard task of determining an asymptotic behavior of a two-point 
correlation function by an easier one, of determining the quantum-critical scaling of 
the quantum fidelity. However, the role of system’s dimensionality has been left as an 
open problem. Our aim in this paper is to fill up this gap, at least partially, by verifying 
the laws of quantum-critical scaling theory of quantum hdelity in a 2D case. To this 
end, we study correlation functions and quantum hdelity of 2D exactly solvable models, 
which are interacting, quasifree, spinfull, lattice-fermion models. The considered 2D 
models exhibit new, as compared with ID ones, features: at a given quantum-critical 
point there exists a multitude of correlation lengths and multiple universal critical 
indices u, since these quantities depend on spatial directions, moreover, the indices 
V may assume larger values. These facts follow from the obtained by us analytical 
asymptotic formulae for two-point correlation functions. In such new circumstances 
we discuss the behavior of quantum hdelity from the perspective of quantum-critical 
scaling theory. In particular, we are interested in hnding out to what extent the 
quantum hdelity approach may be an alternative to the correlation-function approach 
in studies of quantum-critical points beyond ID. 


1 Introduction 

In recent years, quantum phase transitions and quantum-critical phenomena constitute a sub¬ 
ject of great interest and vigorous studies in condensed matter physics. Both, experimental 
and theoretical developments point out to the crucial role that quantum phase transitions 
play in physics of frequently studied high-U superconductors, rare-earth magnetic systems, 
heavy-fermion systems or two-dimensional electrons liquids exhibiting fractional quantum 
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Hall effect [1], [2], Quantum-critical phenomena have been also observed in exotic systems 
as magnetic quasicrystals [3] and in artihcial systems of ultracold atoms in optical lattices [4]. 
The so called classical, thermal phase transitions originate from thermal fluctuations, a com¬ 
petition of internal energy and entropy, and are mathematically manifested as singularities 
in temperature and other thermodynamic parameters of various thermodynamic functions, 
and such characteristics of correlation functions as the correlation length, at nonzero temper¬ 
atures. In contrast, quantum phase transitions originate from purely quantum fluctuations 
and are mathematically manifested as singularities in system parameters of the ground-state 
energy density, which is also the zero-temperature limit of the internal energy density. Nat¬ 
urally, singularities of thermodynamic functions appear only in the thermodynamic limit. 
The importance of quantum phase transitions for physics and the related wide interest in 
such transitions stems from the fact that, while a quantum phase transition is exhibited 
by ground states, hence often termed a zero-temperature phenomenon, its existence in a 
system exerts a great impact on the behavior of that system also at nonzero temperatures. 
A quantum-critical point gives rise to the so called quantum-critical region, which extends 
at nonzero temperatures, in some cases up to unexpectedly high temperatures [5], [6]. 

Theoretically, quantum phase transitions can be studied in quite complex quantum sys¬ 
tems by qualitative and approximate methods, or in relatively simple but exactly solvable 
models by means of analytic methods and high-accuracy numerical calculations [1]. Nat¬ 
urally, for the purpose of testing and illustrating general or new ideas the second route 
is most suitable. Traditionally, this route involves studying the eigenvalue problem of a 
Hamiltonian, the ground state and excitation gaps, determining quantum-critical points and 
symmetries, constructing local-order parameters, calculating two-point correlation functions 
and their asymptotic behavior at large distances and in vicinities of quantum-critical points, 
with correlation lengths and the universal critical indices v that characterize the divergence 
of correlation lengths on approaching critical points. Carrying out such a programme is 
a hard task, which has been accomplished only in a few one-dimensional models. Among 
those models, there are quantum spin chains as the isotropic and anisotropic XY models in 
an external transverse magnetic held, including their extremely anisotropic version-the Ising 
model [1]. Only in one dimension those models are equivalent to lattice gases of spinless 
fermions, which can exactly be diagonalized, and exact results concerning the phase diagram, 
quantum-critical points, correlation functions and dynamics have been obtained (concerning 
XY model see [7], [8], [9], concerning the Ising model see [10], and for both models [11]). 
Needless to say that parallel results for a higher-dimensional model are desirable; this is the 
hrst motivation of our investigations presented in this paper. 

In the last decade, fresh ideas coming from quantum-information science entered the held 
of quantum phase transitions. One of them is the so called quantum-hdelity method. The 
carried out so far studies of particular models show that using this method it is possible to 
locate critical points [12], [13], [14] and to determine the correlation lengths and associated 
universal critical indices v [15], [16], [17]. This is achieved by studying (typically by numerical 
methods) quantum-critical scaling properties, with respect to the size of the system and the 
parameters of the underlying Hamiltonian, of the quantum hdelity of two ground states in 
a vicinity of a quantum-critical point. The index v is extracted from numerically obtained 
plots of hdelity via the quantum-critical scaling laws of quantum hdelity, which have been 
derived by renormalization group arguments [18], [19], [20], [15]. The most comprehensive 
results concerning the hdelity approach have been obtained for one-dimensional quantum 
spin systems in a perpendicular magnetic held [15] (the case of Ising model), [16] (the case 
of XY model), [17] (the case of a quasifree, pairing, lattice-fermion model). These results are 
very promising: except a vicinity of a multicritical point, the hdelity approach works hne. 
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Since the task of determining quantum-critical scaling properties of the quantum hdelity is 
definitely much easier than the task of calculating large-distance asymptotic behavior of two- 
point correlation functions, one is tempted to consider the hdelity method as a substitute for 
the standard correlation-function approach. All that we have said above makes it desirable to 
verify the laws of quantum-critical scaling of hdelity and the ehectiveness of hdelity approach 
in dimensions higher than one, where new features, not encountered in one-dimensional 
models may appear; this is the second motivation of our investigations reported in this 
paper. For that task we need an at least two-dimensional exactly solvable model, whose 
ground states, quantum-critical points with critical indices in their vicinities, correlation 
lengths, and analytic expressions for hdelity can be determined. 

To go beyond the one-dimensional case, we consider lattice-fermion models which orig¬ 
inate from the two-dimensional model of d-wave superconductivity proposed by Sachdev 
[21] (see also [1]), which are spinful pairing models represented by quadratic Hamiltonians. 
In many cases of physical interest the quadratic Hamiltonians are obtained by means of a 
mean-held approximation applied to quartic Hamiltonians of systems with two-body interac¬ 
tions and their parameters are related by self-consistency equations. While our Hamiltonians 
are also quadratic, their parameters are independent—not related by mean-held equations. 
General, mathematical considerations of some classes of such models, but without specifying 
hopping intensities or coupling constants, which therefore do not reach such subtleties as 
quantum-critical points or critical behavior of correlation functions, can be found in [22], [23]. 
For translation-invariant hopping intensities and coupling constants the considered models 
are exactly solvable in any dimension, that is, in particular it is possible to derive analytical 
formulae for quantum hdelity and correlation functions of hnite systems and then in the 
thermodynamic limit, where boundary conditions play no role. To limit further the great 
variety of possible models, we restrict the hopping intensities to nearest neighbors while the 
dimensionality is set to D=2. The underlying lattice is chosen as a square one with hop¬ 
ping intensities invariant under rotations by 7r/2, while the interactions of our systems are 
required not to extend beyond nearest neighbors and to be either invariant under rotations 
by 7r/2 (the symmetric model) or to change sign after such a rotation (the antisymmetric 
model). 

It is worth to mention here another class of models, which might be of interest in the 
context of this paper, known as reduced BCS models of superconductivity and superfluidity, 
see for instance [24] and references quoted there. While those models are exactly integrable, 
analytical formulae for quantities of interest, such as correlation functions, are not avail¬ 
able, except at the thermodynamic limit, where they coincide with the results of mean-field 
theories. For finite systems, all the quantities we are interested in, ground-states, quantum 
fidelity, ground-state correlation functions, are given in terms of only numerically available 
solutions of a set of coupled, nonlinear, algebraic equations, whose number amounts to the 
number of degrees of freedom. Therefore, despite some interesting features, for instance the 
ground state is not a BCS-like state, those models are not suitable for the kind of studies 
reported in this paper. 

The general plan of the paper is as follows. In section 2 we define the two models stud¬ 
ied in this paper, the symmetric model and the antisymmetric one, and give closed-form 
formulae for two basic ground-state two-point correlation functions, which are then used to 
dehne order parameters of those models. Later on, when specific asymptotic behavior of 
two-point functions is discussed, only the gauge invariant function, that is the offdiagonal 
matrix element of the ground-state one-body reduced density operator is taken into account. 
The purpose of the next section 3 is to present the quantum fidelity method of investigating 
quantum-critical points. In particular we provide the known quantum-critical scaling laws 
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obeyed by fidelity in a vicinity of a critical point, and express the fidelity of the models 
considered in the paper as a Riemann sum of an analytically given function of two variables 
- the components of quasimomentum. Then, in sections 4 and 5 we carry out our program 
of confronting predictions of quantum-critical scaling theory of fidelity with exact results 
obtained for 2D pairing models, the symmetric and antisymmetric ones. This program con¬ 
sists of two stages. Since a critical scaling theory is concerned with correlation lengths and 
critical indices, in the first stage analytic results for the spatial direction-dependent behavior 
of the gauge-invariant two-point correlation function at sufficiently large spatial distances 
and sufficiently close to the quantum-critical points exhibited by our models are highly desir¬ 
able. Having such results, one is able to directly infer exact expressions, in terms of system’s 
parameters, for spatial direction-dependent correlation lengths and exact values of the cor¬ 
responding critical indices. We provide such results in the paper; they are excerpts from 
our article [26], where comprehensive studies of the gauge-invariant two-point correlation 
function of our models have been carried out. In the second stage, we are concerned with 
the variations of quantum fidelity against the system’s linear size or against the distance to 
a critical point, sufficiently close to a quantum-critical point. Having expressed the fidelity 
as a Riemann sum, we generate suitable high-accuracy numerical plots that reveal those 
variations. Then, our discussion concentrates on answering two questions. First, can we 
read off from those plots the values of the spatial direction-dependent correlation lengths? 
Second, can we identify regions where the behavior of fidelity matches the known critical 
scaling laws with the known exact values of related critical indices? Anticipating the results 
to be presented, we say only that the answer depends on the kind of a critical point and the 
values of the associated critical indices v. As compared to the results of so far carried out 
analogous studies of ID cases [15],[16],[17], the answer is more likely to be negative. Finally, 
in section 6, we summarize our results and draw conclusions. 

2 The models, their ground-states and ground-state 
correlation functions 

We consider a D-dimensional spinful fermion model, given by the Hamiltonian, 

Li^a 

where a\^, stand for creation and annihilation operators, respectively, of a spin 1/2 
fermion, whose spin projection on a quantization axis is a= ± 1 in units of h/2, in a state 
localized at site l={li, ■ ■ ■ ,Id) of a D-dimensional hypercubic lattice. The edge of the lat¬ 
tice in the direction given by the unit vector e^, i=l ,..., D, whose m-th component is 
consists of Li equidistant sites, labeled by li=0, 1,..., T —1, lj=0 for j^i. In all the consid¬ 
erations that refer to finite systems, special boundary conditions, specified below, are chosen. 
The sums over l,i in (1) amount to the sum over pairs of nearest neighbors, with each pair 
counted once. The real and positive parameter t is the nearest-neighbor hopping intensity, 
p-the chemical potential, J-the coupling constant of the gauge-symmetry breaking interac¬ 
tion, and Aj, i=l ,... ,D, stand for direction-dependent, in general complex, dimensionless 
parameters. In distinction to [21], [1], where Hamiltonian (1) together with equations relat¬ 
ing Aj with /i and J was studied in a context of a mean-field approximation to a kind of t-J 
model, here Aj are free parameters independent of /i and J. Naturally, we can express the 
parameters fi and J in units of t, while the lengths of the underlying lattice in units of the 
lattice constant, preserving the original notation. 


+ h.c.) - 


L t 


J 


D 




aAioja 


t J 

l,cr^l-l-ei,—cr 


+ h.c. 


( 1 ) 
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We note that Hamiltonian (1) is not gauge invariant unless J=0. It is also not hole- 
particle invariant unless /i=0 and A*, i=l,..., D, are real. The latter condition can be 
assumed to hold without any loss of generality, since Hamiltonian (1) with any complex A 
is unitarily equivalent to that with A replaced by |A|. 

Imposing, independently in each direction e*, i=l, ... ,D, periodic or antiperiodic bound¬ 
ary conditions, Hamiltonian (1) can be simplified by passing from the site-localized to the 
plane-wave basis labeled by suitable wave vectors (quasimomenta) k whose components are 
denoted ki, i=l ,..., kf), 

H = -J^COS h + h.c.) , (2) 

k,i 

where Ek stands for the dispersion relation of the hopping term, 

Ek = '^coski - fi. (3) 

i 

In what follows we set Li=L, with even L. Then, the components of wave vectors k can be 
chosen to assume the values ki=7i{2li—Li+2) /Li in the case of periodic boundary conditions 
and the values ki=TT{2li—Li+l)/Li in the case of antiperiodic boundary conditions. Such 
Hamiltonians can readily be diagonalized by means of the Bogoliubov transformation. The 
dispersion relation of quasi-particles reads 

Ek + {^k — Ek ), ( 4 ) 

k 

where J2k ~ is the ground-state energy, and Ek, given by 


Ek — 



i cos ki 


(5) 


are the single quasi-particle energies. For a suitable choice of boundary conditions specified 
above, as long as our system is finite the excitation energies Ek remain strictly positive: 
F'fc > 0 for all values of k, and this is assumed to hold in the sequel. Specifically, to 
avoid closing down of the gap in hnite systems, represented by the two-dimensional models 
discussed in the sections that follow, we impose periodic boundary conditions along one axis 
and antiperiodic ones along the orthogonal axis. 

The Hamiltonian (1) preserves parity; therefore without any loss of generality we can 
restrict the state-space to the subspace of even number of fermions. In this subspace, the 
state |0)gp - the quasi-particle vacuum of an unspecified (but even) number of fermions, 
defined by 

|0)<?P = + '«^fc4^^c4,|)|0), (6) 

k 

where |0) is the fermion vacuum, with Uk real and positive. 



and, in general, complex Vk, 


l^fcl = 



^k 

Ek 


arg Vk = arg 




cos ki 


(7) 

( 8 ) 
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is the eigenstate of ( 2 ) to the lowest eigenenergy, (sk — Ek). As long as -Efc > 0 for all 
values of fc, the unique ground state | 0 )gp is the vacuum of elementary excitations (quasi¬ 
particles). However, on passing to the thermodynamic limit, when the system’s linear sizes 
in all directions tend to inhnity, the minimum of Ek over k (i.e.the excitation gap in the 
spectrum of quasi-particles) may approach zero at special values of the chemical potential 
H and the coupling constant J. Those special points in the (p., J)-plane, where Ek, as a. 
function of continuous wave vector k, vanishes, constitute the quantum-critical points at 
which the system undergoes continuous quantum phase transitions. 

All the correlation functions of considered systems can be expressed in terms of two 
basic two-point correlation functions. Since we are interested only in ground-state correlation 
functions, taking into account the lattice-translation invariance of our system these two basic 
two-point correlation functions can be chosen as follows: 

5 p( 0 1 Oq |„-Or ,(T 10)gp and gp(0|cio,(T®r,—crlO)^^, (9) 

with some a. The hrst correlation function, gp(0|aQ „ar,a\0)qp, is gauge and spin-flip invariant; 
for r 7 ^ 0 it represents offdiagonal matrix elements of the ground-state one-body reduced 
density operator. Expressing aj,^, by the creation and annihilation operators of quasi¬ 
particles (the Bogoliubov transformation) we get 

gp{0\al^ar,a\0)qp = \vk\‘^ exp tkv = exp ifcr, ( 10 ) 

k k ^ 

which, upon using the invariance of Sk and Ek with respect to reflections of k in coordinate 
axes, in the thermodynamic limit becomes 

lim qp{Q\al ar,a\^)qp = G{r) = / dk-^ TT cos kjrj . ( 11 ) 

L^oo 271^ Lk y 


Choosing the spin projection a = -|-1, the second correlation function, measuring the degree 
of gauge-symmetry breaking, amounts to 


qp( 0 |ao,+ar-| 0 )gp = -^y^MfcUfcexp -ikr 


2L^ 


E 


J Ai cos ki 


E, 


k 


exp ikr, ( 12 ) 


k k 

which, by the above arguments, in the thermodynamic limit becomes 


lim qp{0\ao,+ar-\0)qp = h{r) = - 7 ^ / dk 

L -,00 2tI^ J[0,7r]O 


J Yji A* cos k. 
Ek 


D 


- COS k 

i=i 


Aj- 


(13) 


Both the above dehned two-point correlation functions are used to dehne the order pa¬ 
rameters in the ground-state phase diagrams of the models considered below. The two-point 
correlation functions decay with increasing distance |r| between the points. In gapped phases 
their decay is dominated by an exponential factor, exp(—|r|/^), which dehnes the correlation 
length f. If additionally (p., J)-points approach a quantum-critical point, i.e. the distance 6 
between them tends to zero, then f diverges as 6~^, which in turn dehnes a universal critical 
index u associated with a particular quantum-critical point. For D>1, two-point correlation 
functions, and hence correlation lengths f as well as the critical indices u may depend on 
spatial direction of vector r . 
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In what follows we shall study only the case of D=2. Clearly, all the above D-dimensional 
expressions can be adapted to the 2D case by setting Ai=ki=ri=0 for i > 2. Note however 
that even with the restriction D=2, due to the freedom in choosing the relation between 
the parameters Ai and A 2 , formula (1) represents a great variety of models. In this paper 
we limit our considerations to two cases only. First we set Ai=A 2 =A, what results in 
the interaction term invariant under rotations by 7r/2; the corresponding model is dubbed 
symmetric. Then, we choose Ai= — A 2 =A, what results in the interaction term changing 
sign under rotations by 7r/2; the corresponding model is dubbed antisymmetric. In all the 
considerations below, that refer to the symmetric or the antisymmetric model, we make 
the identification J|A|=J. In both the cases the correlation functions of our systems are 
invariant not only with respect to lattice translations but also with respect to rotations by 
7r/2. 

In order to determine the correlation length ^ and the critical index in a vicinity of 
some quantum-critical point we study the so called doubly-asymptotic behavior of the gauge- 
invariant correlation function G{r), defined in(ll); the provided in this article analytical 
results are excerpts from more comprehensive studies presented in [26]. Clearly, for a fixed 
spatial direction and the parameters Aj, G{r) depends on three parameters: |r| - the distance 
between the two points of the correlation function, the chemical potential fi and the coupling 
constant J. The so called doubly-asymptotic region is characterized as follows: for a fixed 
but sufficiently large |r| the (/r, J)-points approach a quantum-critical point along a specihc 
path. The choice of those paths depends on the kind of considered critical points. The critical 
points of the symmetric model that are located at the J-axis, with | J|>0, and those of the 
antisymmetric model that are located on the lines parallel to the J-axis are approached along 
/i-paths that are parallel to the /i-axis. Then, the critical points of the symmetric model 
that are located at the /i-axis are approached along J-paths that are parallel to the J-axis. 
Finally, the multicritical point of the symmetric model is approached along the 45°-path. 

As mentioned above, a novel feature of two-dimensional models, in comparison with the 
one-dimensional case, is that the two-point correlation function G{r) depends not only on 
the distance |r| but also on the direction of r. Expressing r by its Cartesian coordinates, 
r=(ri,r 2 ), we can parameterize directions by the ratio ri/r 2 =n. Then, for a given critical 
point, we can expect n-dependent doubly-asymptotic behaviors of correlations. Unfortu¬ 
nately, the analytic asymptotic formulae for G{r) in offdiagonal directions (i.e. n^l), which 
we have been able to obtain, apply only to points r such that ?7,>?7,o>l or, by symmetry, 
n<nQ^<l, that is for offdiagonal directions which form a sufficiently small angle with the 
axial directions. Therefore, analytic asymptotic formulae in the diagonal direction (i.e. n=l) 
have been derived separately [26]. These formulae define u-dependent correlation lengths 
^offdiag offdiagonal directions satisfying the conditions specihed above and the correlations 

length in the diagonal direction, where the superscript plus refers to the symmetric 
model and minus - to the antisymmetric one. For directions close to the diagonal one, that 
do not satisfy the above conditions, the correlation lengths and related critical indices have 
been determined numerically. Interestingly, our analytical and numerical results show that, 
for each critical point of the symmetric or the antisymmetric model, there are only two kinds 
of universal critical indices u: VoEdia.^ offdiagonal directions and i^diag for the diagonal 

direction. 
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3 The ground-state fidelity and quantum-critical scal¬ 
ing laws 

Let A be a vector whose components are those parameters of the considered system’s Hamil¬ 
tonian that drive a quantum phase transition, and e - a unit vector in the space of those 
parameters. Then, on varying parameter 6 the vectors A -|- 5e scan a neighborhood of A 
along direction e. For given A and 6, the ground-state hdelity at A in direction e, Te(A, 5), 
is the absolute value of the overlap of the ground states |A ± 6e) at the points A ± 6e, 

Te(A,5) = |(A-(5e|A + 5e)|. (14) 

A list of general, system independent, properties of Te(A, 5) can be found in [25]. The 
transition point of a continuous quantum phase transition, i.e. a quantum-critical point, 
denoted Ac, is characterized by the power-law divergence of the correlation length .^(A), as 
the quantum-critical point is approached: ~ |A — Ac|“'^, with u being one of universal 

characteristics of a critical point. Alternatively, Ac can be dehned as the point where the 
gap between the ground-state energy and the energy of the lowest excited state vanishes. In 
reference [12] it was demonstrated that the quantum-critical point Ac can be identihed as 
the value of A at which the minimum of hdelity is realized. However, the hdelity approach 
seeks an answer to a more general question: does the behavior of quantum hdelity in a 
neighborhood of a quantum-critical point encode not only the location of that point but 
also some universal properties of the underlying quantum phase transition? First results, 
pointing towards a positive answer to the raised question by providing some hnite-size critical 
scaling of hdelity, have been obtained by Venuti and Zanardi [18]. 

Scaling theories of criticality are based on the assumption that the behaviour of a system 
is governed by a unique length—the correlation length. In particular, according to hnite-size 
scaling theories, the properties of a system are close to those at the thermodynamic limit, 
if the linear size of the system, L, is much greater than the correlation length .^(A). In 
literature such a system (or regime) is called the thermodynamic limit system (regime) or 
oh-critical system (regime). In the opposite case one speaks of a hnite system (regime) or 
a quasicritical system (regime). However, discussing hnite-size scaling properties of ground- 
state hdelity, Te(A, 5), we deal with two ground states of a system. Then, one of them 
may correspond to a quasicritical or critical regime while the other to an oh-critical one. 
In the context of the quantum-critical-scaling theory of hdelity, it is expected that the role 
of the unique characteristic length of the system, that diherentiates between a quasicritical 
system and an oh-critical one, is played by the ehective correlation length .^e(A, 6), which is 
dehned as the smaller of the two correlation lengths .^(A ± 6e) [16]. To turn attention to 
this fact and to emphasize that in all the considerations below the systems are hnite, we 
introduce a somewhat nonstandard terminology. If we call the system macro¬ 

scopic; in the opposite case the system is called small. The crossover between small-system 
and macroscopic-system regimes occurs, when the linear size of a system, L, satishes the 
crossover condition: 

L/U\5)r^l. (15) 

There are numerous papers devoted to critical scaling of small-system hdelity, see [25], 
[20], and references quoted there. Typically, small-system hdelity can be Taylor-expanded 
in 6, 

y.(A.i) = l-h^(A)+..., 
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( 16 ) 


where the first order term vanishes because of the symmetry of fidelity in 6 at zero. The 
coefficient of the second order term, Xe(A), is known as the fidelity susceptibility. One expects 
some universal scaling properties of Xe(A), provided A is sufficiently close to a quantum- 
critical point Ac, where the correlation length diverges: .^(Ac ± Se) ~ Fairly general, 

model-independent, arguments provide us with finite-size scaling of the fidelity susceptibility 
at Ac [191, 201: 

XeiK) ~ L^'\ (17) 

or equivalently, in the small-system regime 

-lnTe(Ac,(5) ~ (18) 

Let us note here that in vicinities of some quantum-critical points fidelity oscillates on varying 
L, with an amplitude that is particularly large, close to one, in the small-system regime [16], 
[17]. In such cases, the hdelity susceptibility is not well defined. However, the small-system 
scaling law (18) may still hold but in a generalized sense [17]. Specifically, it is the envelope 
of the minima of — lnTe(Ac, (5) that scales according to (18). 

In the macroscopic-system regime, quantum phase transitions have been studied by means 
of the so called hdelity per site, a quantity whose logarithm is equal to N~^ In Te(A, 5), where 
N = is the number of sites in a H-dimensional system [13], [14]. However, critical scaling 
of a macroscopic-system hdelity has been considered only very recently by Rams and Damski 
[15],[16]. These authors have found that, while for small-systems the hdelity scaling is totaly 
insensitive to the way the critical point Ac is approached by the points A ± (5e (i.e. for 
instance, whether they are located on one side of the critical point or on the opposite sides), 
in the case of macroscopic-system the way of approaching the critical point matters. To 
make this explicit. Rams and Damski substituted Xc + c6e for A. By choosing the value of 
the parameter c, the above mentioned location of the two points Ac -|- c6e ± 6e with respect 
to the critical point can be controlled. If |c|>l or |c|<l, then both points are located on one 
side or on opposite sides of Ac, respectively. If |c|=l, then one of the points coincides with 
Ac. Now, the above mentioned independence of the small-system-fidelity scaling on the way 
the critical point Ac is approached by the points A ± (5e can be expressed as follows: 

- In Te(Ac + cSe, 6) ~ (19) 

In contrast to the small-system case, the fidelity scaling law for macroscopic systems, derived 
by Rams and Damski [15], makes the dependence on parameter c explicit. Provided that 
the thermodynamic limit of In Te(A, 5) does exist, it reads 

-lnTe(Ac + c5e,(5) ~ \6\^’^NAe{c), (20) 

where Af>{c) is the scaling function. 

It should be emphasized that the small-system scaling law (19) as well as the macroscopic- 
system scaling law (20) have been derived, using critical-scaling theory arguments, under two 
conditions. The first one is that there is only one characteristic length scale in the underlying 
system, which discriminates between small systems and macroscopic systems. This charac¬ 
teristic length is identified with the effective correlation length ^e(Ac + c5e, 5)=^e(Ac, 5). The 
second one is that the strict inequality Dz/< 2 holds true [19], [20], [15]. It is worth to mention 
that if this condition is satisfied, then in the small-system regime Xe(Ac) or — In Te(Ac, 5), 
formulae (17), (18), respectively, scale with system’s linear size in a superextensive way. 

Let |0)qp and |0)qp be two ground states, the first one for pairs (/i, J), and the functions 
i7fc, the second one for pairs (/I, J), and the functions dfc, E^. As a result of the product 
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structure of the ground states, the quantum hdelity for these states has also a product 
structure, 

lijp ( 0 | 0 ) 

9pl n I {ukUk + \vk\\vk\expi{aTgVk - argu*,)) |. (21) 

k 

After using (7), ( 8 ) the hdelity of two ground states ( 6 ) assumes the form 

l.(0|6)„| = n W. m = 5 (i + ■ P2) 

where the sum of cosine functions has to be taken in the case of the symmetric model and 
the difference—in the case of the antisymmetric one. 

Let us adapt the general notation introduced in the begining of this section to the con¬ 
sidered models. As the location of critical points is uniquely determined by pairs (p, J), 
we set A=(p, J), hence |A) = |0)gp. Then, in formula (22) for hdelity, the functions Ek, 
given by (3) and (5), respectively, are calculated at A^ -|- (c — l)6e, while ik and Ek - at 
Ac + (c-|- l)Se. Finally, we set | 5 p( 0 | 0 )gp|=Te(Ac, 5). 

Considering quantum-critical scaling of hdelity, we shall study numerically the sum 

-i5^1n/(fc) = -lnTe(Ac,<5), (23) 

k 

as a function of parameter 6 for hxed system size N, or vice versa, in neighborhoods of 
various critical points, in small- and macroscopic-system regimes. We recall that the values 
of ki are obtained with periodic boundary conditions and the values of /c 2 with antiperiodic 
ones or vice versa. In all the considered cases the function In f{k) is either continuous in 
the whole square [ 0 , 7 r]^ or it has an integrable singularity at some k (a discontinuity or 
a logarithmic divergence). Therefore, in all the considered cases the limit iV—)-oo of the 
Riemann sum corresponding to (23) does exist, 

lim -7VMnTe(Ac,(5) = 7 ^ f dfcidfcs (-In/(fci, ^ 2 )). (24) 

Consequently, for given sufficiently small 5 and sufficiently large N 

- lnTe(Ac,(5) ^ 7 ^ / d/cid/c 2 (-ln/(/ci,/c 2 )) (25) 

2[0,7r]2 

approximately, that is in a macroscopic-system regime — lnTe(Ac,(5) scales with the system 
size as N. 

Any study of critical scaling involves specifying a critical region, that is a critical point 
and its neighborhood. The quantum-critical points considered in our paper are displayed in 
Figs. 1 and 7. As for neighborhoods, we have chosen line neighborhoods, each one specihed 
by a unit vector e and a range of parameter 5, which are scanned by vectors Ac -|- (c — l)6e 
and Ac -|- (c + l)6e on varying 6. Without any loss of generality only (5>0 is considered. 

Our aim in the sequel is to confront the predictions of quantum-critical scaling theory 
for hdelity with exact results, to hnd out limitations and advantages of hdelity approach in 
investigations of quantum-critical points in dimensions D>1. To the best of our knowledge, 
this task has never been carried out. We note that the case D=1 has already been extensively 
studied in [15], [16], [17] and the results are promising. 

We would like to learn to what extent the quantum hdelity is useful for determining 
the correlation lengths and related universal critical exponents u in D=2 systems, that is 
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in cases where there exists a multitude of those quantities of interest, due to their spatial- 
direction dependence, and the inequality Du<2 may not be satished. Thus, in such cases 
one or both the basic assumptions of the quantum-critical scaling theory of hdelity are 
violated. For this purpose, we consider different kinds of critical points exhibited by the 
considered systems and calculate — lnTe(Ac,5) as a function of 6, keeping the linear size L 
hxed or vice versa. Then, we make an attempt to determine the intervals of 6 or L, where 
— lnTe(Ac, S) obeys a power law. After that, we try to identify the regime of small system, 
the one of macroscopic system, and the characteristic length, or lengths, that discriminates 
between the small- and macroscopic-system regimes - the location and extent of the crossover 
regime. In connection with this issue we note here that according to our results concerning 
correlation lengths, summarized in sections that follow, for a given linear size of the system, 
the effective correlation length increases, we move towards the regime of small system, by 
decreasing sufficiently |5|. In the opposite case the effective correlation length decreases and 
we move towards the regime of macroscopic system. Finally, in regions where — lnTe(Ac, <5) 
obeys a power law, using the scaling laws (19) and (20) we read off from the plots of hdelity 
the values of u and compare them with the obtained exact values given in the displayed 
phase diagrams. 

While the two-point correlation function G{r), in particular its large-distance asymptotic 
behavior, and consequently the dehned above effective correlation length ^e(Ac, 5) depends 
on spatial directions, hdelity Te(Ac, S) does not even ’’know” what a spatial direction is. This 
fact rises the interesting question of the relation between the mentioned above characteris¬ 
tic length, or lengths, and a multitude of spatial-direction-dependent, ehective correlation 
lengths .Ce(Ac, 5). Anticipating our discussion of results presented in the sections that follow, 
we say only that we have demonstrated that the mentioned multitude of ehective correlation 
lengths is rehected in the existence of more than one crossover regimes whose location is 
given by the ehective correlation lengths in specihc spatial directions. 


4 The case of symmetric model 

We can distinguish four ground-state phases labeled by two order parameters, Oi and O 2 , 
dehned as 

Oi = G(0)-^, O2 = -A*h{l,0). (26) 

To locate the quantum-critical points we solve the equation i?fc=0. An inspection of formula 
(5) leads immediately to the conclusion that there exist wave vectors k such that Ek=0 if and 
only if /i=0 and J is arbitrary (the whole J-axis in the (/i, J)-plane) or J=0 and /iG[—2,2] 
(the closed interval [—2,2] of the /i-axis). There are two critical end points (±2,0) and a 
multicritical point (0, 0). The ground-state phase diagram of the symmetric two-dimensional 
system is shown in Fig. 1. In all the analytic asymptotic formulae presented below, a vicinity 
of the multicritical point (0, 0) is excluded. In particular, for /r —)■ 0 the J-coordinates of 
/i-paths have to be away from zero; analogous condition applies to J-paths. The case of 
multicritical point will be discussed separately. 

In the stripe |/i|<2 of the (/i, J)-pIane, but excluding the /i=0 and J=0 lines, the large- 
distance asymptotic behavior of G{r) is: 

G(r', r') ^ -sgn(/i)^ ^ cos(0r' ± 0), (27) 
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Figure 1: Phase diagram of the symmetric two-dimensional system in the (/i, J)-plane. The 
set of quantum-critical points consists of the J-axis and the closed interval [—2, 2] of /i-axis 
- thick lines. Those lines constitute also phase boundaries of the four phases, labeled by the 
order parameters Oi and 02 - Double arrows indicate the types of neighborhoods of critical 
points, in which the asymptotic behaviors of G(r) are studied, except neighborhoods of the 
multicritical point. The universal critical indices ly in those neighborhoods, whose values are 
given by the arrows, depend in general on spatial direction, whether it is diagonal (i^diag) or 
offdiagonal (z/oMiag)- 


in the diagonal direction, and 


G(ri,r2) - 


Cr f fi 


2J2 y/Vl + exp (-ry^Lg) 


271 \1 + J^ 




COs(r6'ofTdiag + 0), 


with 


r = ri-I- Cr = 

1 


1, if /i>0, 

(_l)(n+r 2 +l)^ if/i<0. 


d(+) 

^offdiag 


n 


1/2 


and. 


e, 


offdiag 


1 + m? 

2 \ 1/2 
n ' 


1 -|- 


1 

6 ^ 


^ -^^2 J , 


(28) 


(29) 

(30) 


(31) 


provided that the points (ri,r 2 ) become remote from the origin along a ray ri/r 2 =n=const 
[26]. We recall that by symmetry the large-distance asymptotic behavior of G{r) in offdiag¬ 
onal directions is the same for ?7,>?7,o>l and for n<nQ^<l (specihcally we found that one can 
choose no=3). However, the above formulae that refer to offdiagonal directions hold only for 
n>nQ>l. The formulae (27) and (28) dehne the diagonal, and offdiagonal, 

■^ofTdiag’ correlation lengths, respectively. The parameters 9, .^i, ,^ 2 , 9i and 62 will be 
expressed by p and J in vicinities of critical points in subsections that follow, hence the 
behavior of G{r) in doubly-asymptotic regions will be specihed. Now, we are ready to carry 
out our programme for specihc types of critical points exhibited by the symmetric model. 
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Figure 2: (Color online) The symmetric model, critical point /i=0 and J=l, 6=fi. Plots of 
— InT(i,o)((0,1),(5) for three values of c: c=0 - red line, c=l - green line, c=2 - blue line. 
Left panel: plots of — lnT(i,o)((0,1), 10“^) versus L. Right panel: plots of — lnT(i_o)((0,1), S) 
versus 6 for L=10"^. Both plots are in doubly logarithmic scale. Here and in all the hgures 
below, black dashed-dotted straight lines indicate the power-law scaling. The uncertainties 
in their slopes have been obtained by the least-square htting. 


4.1 Critical points at the line /u=0 

In this subsection we consider (ju, J)-points approaching along a p-path a point belonging 
to any one of the two half lines of quantum-critical points, given by ju=0 and |J|>Jo>0, 
for some Jq. Then, in terms of /i and J the parameters determining the large-distance 
asymptotic behavior of G{r) in (27) and (28) are given by (see [26]) 


1 _ \^^■A 

^ 1 + J2 ’ 


6 ^ TT — 


l/^l 

1 + J2’ 


(32) 


1 

c 


2|/i| 


vm^ 


sm 


arctan | J| 



6i ^ TT — 26*2 


(33) 


^2 ~ y 2 V 1 + arctan I J|^ . (34) 

In the left and right panels of Fig. 2 we show two typical plots of — In T(i o) ((0,1), 5); in the 
left panel hdelity is plotted versus the system’s linear size L while in the right one - versus 
the deviation from the critical point 6=fi. Three regions of different behavior of hdelity, 
separated by two crossover regions are well visible. Our hrst goal is to relate the locations 
of the crossover regions with the correlation lengths of the system. It appears that in the 
left panel plot of hdelity, made for hxed 5=10“^, on increasing L the hrst crossover occurs in 
vicinity of the ehective correlation length 'C(i,o)((0,1), S) in the axial direction (later denoted 
■CaxiaK*^)), while the second one - in vicinity of the ehective correlation length ^(i,o)((0,1), <5) in 
the diagonal direction (denoted ^diag(<5))- With the two ehective correlation lengths, ^axiai(<5) 
and ^diag(^), we can associate two ehective deviations S, i.e. (5axiai(L) and (5diag(T), dehned 
as the solutions of the equations ^axiai(<^)=-f" and ^diag(<5)=T, respectively. It appears that 
in the right panel plot of hdelity, made for hxed L=10^, on increasing 6 the hrst crossover 
occurs in vicinity of the ehective deviation ^axiai(-f") while the second one - in vicinity of 
the ehective deviation 6diag{L). Having identihed the locations of the crossover regions with 
ehective correlation lengths in specihc directions, we would like to assign a meaning for them, 
dehne them using the set oh all the spatial direction-dependent correlation lengths. To this 
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end, let us denote by ^n(<5) the effective correlation length ^(i,o)((0,1),^) in spatial direction 
given by n>l and let n 2 >ni>l. From formulae (30), (32), and (33) (see also Figs. 15, 16 
and comments in section 5 in [26]), one easily infers that for sufficiently small 5 the following 
inequalities hold true: 

eaxial(5) < < fdiag(5), (35) 

provided ?7,i>no>l; moreover, ^axiai(5)= hm„^oo'Cn(<5). However, our numerical results pre¬ 
sented in [26] support the hypothesis that inequalities (35) hold true for any n>l, and then 
^diag(<5)= lim„_,.i ^n(^)- Let 5n{L) be the effective deviation in direction n, i.e. the solution of 
the equation ^^{5)=L for some sufficiently large L. Since ^n(^) is a decreasing function of <5, 
inequalities (35) between effective correlation lengths in different directions imply analogous 
inequalities between the corresponding effective deviations: 

^axial(L) < 5^12 (-L) < ^ni (L) < 5diag(L), (36) 

for sufficiently large L, where 5axiai(L)= hm„_,.oo (5„(L), (5diag(L)= hm„_,.i 6 n{L), and vice versa. 
From (35) one concludes that for sufficiently small deviation 6 

^axial(*^) lllf ^n((5) = ^lower(*5) aud C^diag(*^) SUp^[i(*5) = C^upper(*5)) (^'^) 

n>l 

where we defined the lower effective correlation length Slower (^) and the upper effective cor¬ 
relation length Supper(^)- Analogously, for sufficiently large L, (36) implies that 

*^axial(L) inf Sn(L) = (5lower(L) aud <5diag(L) SUp(5fi(L) = (5upper(L), (38) 

n>l n>l 

where we defined the lower effective deviation (5iower(L) and the upper effective deviation 

*^upper (L). 

Summing up, our numerical results, in particular those displayed in Fig. 2, show that 
there are two characteristic lengths, .^iower(5) and Supper (<5), or equivalently two characteristic 
deviations from the critical point, Slower(L) and Supper(L), that mark the crossover regions in 
the behavior of — ln5'(i_o)((0,1),^) versus L or 5, respectively. Remarkably, in spite of the 
fact that the hdelity does not depend explicitly on spatial directions n, a finger-print of the 
dependence of correlation properties on spatial directions can be seen in the behavior of the 
fidelity as a function of L or 5. 

We can naturally identify the range of sufficiently large but smaller than Slower (^) values of 
L, or the range of sufficiently large but smaller than Slower (L) values of <5, as the small-system 
regime. Then, the ranges of L and <5, given by double inequalities flower(5) <L<.^upper(<5) and 
(Slower (-L)<5<Supper (L) cau be identified as the mesoscopic-system regimes. Finally, the ranges 
of L and S above Supper (<5) or Supper (L), respectively - as the macroscopic-system regimes. 

Having defined the characteristic lengths, .Ciower(<5) and ^upper(<5)5 and characteristic devia¬ 
tions from the critical point, Slower(L) and Supper(L), which mark the crossover regions in the 
behavior of fidelity, and then the regimes of small- and macroscopic-system, we can describe 
our numerical results, from the perspective of the scaling laws (19), (20). We note that the 
divergence of Slower (<5)=^axiai(<5) as (5—j-O is characterized by zZoffdiag=l/2, which satisfies the 
condition Du<2. Then, in the small-system regime of L<^iower(^) the scaling of fidelity 
is consistent with (19) for z/=l/2, which matches r'offdiag- In the macroscopic-system regime 
of L>.^upper((5) we observe the standard (see section 3) scaling of fidelity. On the other 
hand the scaling of fidelity with respect to 5 does not provide any information about the 
exponent v. In the small-system regime of 5<(5iower(-L) we observe the standard (see section 
3) (5^ scaling of fidelity, while in the macroscopic-system regime of (5>5upper(L) the fidelity 
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Figure 3: (Color online) The symmetric model, critical point /i=l and J=0, 6=J. Plots of 
— lnT(o,i)((l, 0), (5) for three values of c: c=0 - red line, c=l - green line, c=2 - blue line. 
Left panel: plots of — lnT(o,i)((l, 0), 10“®) versus L. Right panel: plots of — lnT(o,i)((l, 0), 6) 
versus 6 for L = 10^. Both plots are in doubly logarithmic scale and the black dashed-dotted 
straight lines indicate the power-law scaling. The variable L in the left panel changes by 
4 up to 10 000, but above 10 000 the formula L —)■ 4[1.05L/4] is used, where [■] denotes 
rounding up to the nearest integer (ceiling function). The latter way of sampling L reveals 
the behaviour for sufficiently large L, where for computational as well as presentation 
reasons, continuing to change L by 4 is no longer feasible. This formula has an advantage 
of evenly spacing values of L in logarithmic scale but fails to capture properly the character 
of the oscillations for sufficiently small L. 


behaves in an anomalous, non-power-law way, which does not allow for estimating u via (20). 
We note however, that Supper(<5), which is used to calculate Supper(L), diverges as <5“^ if )-0 
(z/diag=l); which violates the condition Du<2. 

4.2 Critical points in the intervals J=0, 0<|/i|<2 

In J-path neighborhoods of the critical points at the line segments J=0 and 0<|/i|<2, 


1 2|/iJ| 

~ ^4 - /i2 ’ 


9 K, 2 arccos 


l/^l 


(39) 


1 

6 


\J\ 



6i ^ IT — 2 arcsin 


./M + 1 / l/^l 3 - |/i| 2 

V 2 ^2\/2-|/i|2-|/i| 


(40) 


1 ^J/^l - 1 1 
6 ^ 2 6 ’ 




(41) 


for sufficiently small \J\ [26]. 

In those neighborhoods, the effective correlation lengths 'C ('^)5 functions of 6=J and 
n, share the properties that were used in the previous subsection to define the lower and 
upper effective correlation lengths, flower (<5) and .^upper(<5)- Therefore, we can adopt the same 
definitions of these lengths, the associated effective deviations, Slower (L) and Supper (T), and 
the small- and macroscopic-system regimes. In the special case of /i=l, all the hierarchy of 
correlation lengths and the associated deviations collapses, 6ower(5)~'Cupper(<^) for sufficiently 
small S, and ^lower (-f^) ~(Supper (-f^) for sufficiently large L. The critical exponents characterizing 
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the divergence of ^iower(^) and ^upper(^) coincide, r'offdiag=i^diag=l) and violate the condition 
Du<2. In Fig. 3 we show an example of the behavior of fidelity in neighborhoods of the con¬ 
sidered critical points. In the left panel, the plot of — In 5'(o,i)((l, 0), 10“®) versus L exhibits 
pronounced oscillations, whose amplitude is particularly large for L<,^si and decreases with 
increasing L. We remark that the nature of those oscillations has been revealed in [16]. As 
we mentioned in section 3 (see also [17]) on presenting small-system scaling law, in case of 
oscillating hdelity one should consider scaling in the generalized sense. That is, in such a 
case it is the envelope of the minima of — lnTe(Ac, 5) that is the right quantity whose scal¬ 
ing should be studied. In the left panel, one crossover region separating different power-law 
behaviors, in the generalized sense, marked by some ^si(<5), is visible. However, does not 
match 6ower(<5)~^upper(<5), whosc valucs are a lot larger than ^si(<5). If we naively identify 
the regime of T<^si as the small-system regime, then the observed L^-scaling of hdelity im¬ 
plies, via (19), i^=l/2, which does not match the exact value 1 of this exponent. Moreover, 
L^-scaling of — In Te(Ac, 5), consequently superextensivity, interestingly occurs despite the 
fact that the condition of superextensive behavior, which follows from critical-scaling theory, 
Du<2, is violated. For L>^si, identihed as the macroscopic-system regime we observe the 
standard LAscaling. In the right panel, the plot of — lnT(o,i)((l, 0), 5) versus S exhibits also 
one crossover region, separating different power-law behaviors, marked by some 6si, which is 
the solution of the equation ^si{S)=L. The power law (20) applied for 6 well above (5si gives 
the incorrect result z/=l/2. We note that similar results hold for other values of /r, |/i|7^1. 
The only difference is that 6ower((5)<^upper(5), but still 6ower((5)>6i(<5); thus, 6ower((5) does 
not locate the crossover region correctly. 

4.3 The end critical points J=0, |/i|=2 

If (/i, J)-points approach along the J-path one of the two end critical points, |/i|=2 and J=0, 



e ^ 2 |j|h 2 , 

(42) 

V2i7=|. 

01 ^ = 202, 

(43) 


for sufficiently small J [26]. 

In a J-path vicinity of one of the end critical points J=0 and |/i|=2, setting 6=J, 
we can adopt the same dehnitions of ^iower(<5) and ^upper(<5) of the small-, mesoscopic- and 
macroscopic-system regimes as in subsections 4.1 and 4.2. However, using (42) and (43) one 
verihes easily that, for sufficiently small J=6, ^axiai(<5) ~ ^n(^) ~ ^diag(^), hence ^iower(<5) ~ 
Supper (5). Therefore, there is no mezoscopic-system regime, there is only one characteristic 
length, say .^iower(<5), one corresponding (5iower(-h), that discriminates between the small- and 
macroscopic-system regimes. Indeed, in Fig. 4, in the left panel, where hdelity is plotted 
against L, as well as in the right panel, where hdelity is plotted against 6, only one crossover 
region, marked by .^iower(<5) or (5iower(-h)5 respectively, is visible. The critical exponents char¬ 
acterizing the divergence of .^iower(<^) and ^upper(<5) coincide, i^offdiag=^'diag=l/2, and satisfy 
the condition Du<2. In the small-system regime, L<.^iower(^) (left panel), hdelity scales as 

which via (19) gives the correct value of the exponent i^=l/2, while in the macroscopic- 
system regime, L>^iower(<5), the standard L^-scaling is observed. Then, in the small-system 
regime, (5<(5iower(-h), hdelity exhibits the standard (5^-scahng, while in the macroscopic-system 
regime, 5>5iower(-h), it scales as 5, which via (20) gives again the correct exponent z/=l/2. 
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Figure 4: (Color online) The symmetric model, critical point /i=2 and J=0, 6=J. Plots of 
— lnT(o,i)((2, 0), (5) for three values of c: c=0 - red line, c=l - green line, c=2 - blue line. 
Left panel: plots of — lnT(o,i)((2, 0), 10“®) versus L. Right panel: plots of — lnT(o,i)((2, 0), 6) 
versus 6 for L = 10^. Both plots are in doubly logarithmic scale and the black dashed-dotted 
straight lines indicate the power-law scaling. For more details see the text. 



Figure 5: (Color online) The symmetric model, critical point /i=J=0, 6=fi=J. Plots of 

— lnT(i^i)((0, 0), 5) for three values of c: c=0 - red line, c=l - green line, c=2 - blue line. 
Left panel: plots of — In T(i4)((0, 0), 10“'^) versus L. The calculated effective correlation 
lengths ^axiai(<5) and ^diag(^) are too large to be depicted in the hgure. Right panel: plots 
of — lnT(i^i)((0, 0), 5) versus 6 for L = 10^. The effective deviations ^axiai(L) and (5diag(L) 
are well above Smm- Both plots are in doubly logarithmic scale and the black dashed-dotted 
straight lines indicate the power-law scaling. For more details see the text. 

4.4 The multicritical point fi=0=J 

In this case we cannot provide any analytic results concerning the correlation lengths and 
associated critical indices shown in the phase digram of Fig. 1. These quantities have been 
obtained numerically by studying the behavior of G{r) in doubly asymptotic regions, where 
the multicritical point has been approached along the 45°-path {^=J). To determine i^ofrdiag, 
a few spatial directions, including the axial ones, have been taken into account. Both the 
exponents extracted from the asymptotic behavior of G{r), i.e. r'ofrdiag=3/2 and r'diag=2, 
violate the condition Du<2. To the best of our knowledge, /i=0=J quantum-critical point 
is the first instance of quantum critical point whose associated indices v are so large that 
they satisfy the strict inequality Du>2. 

The plots of hdelity are shown in Fig. 5, against L - in the left panel and against 6=fi=J 

- in the right one. Two crossover regions, marked by ^sm, Cmm, with ^sm<^mm, and Ssm, 
Smm: with Ssm<Smm, respectively, are visible. Numerically calculated effective correlation 
lengths .^axiai(<5) and ^diag(<^) are so much larger than ^rnm that they cannot be displayed in 
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the left panel of Fig. 5. The corresponding (5axiai(-f") and (5diag(-^) are larger than 5mra as 
well but they are visible in the right panel. If we naively, ignoring critical scaling theory, 
identify the region of L<^sm (^<^sm) as small-system regime and that of L>^mm {S>6mm) 
- as macroscopic-system regime, we can describe them from the perspective of the hdelity 
scaling laws. We observe the standard L^-scaling of hdelity for L>^rnm and the standard 
h^-scaling for 5<5sm- For 5>5mm no prediction for u can be made, since hdelity scales with 
h in a non-power-law way. For L<^sm, hdelity scales as L^, which gives, via (19), i>=l - the 
value well below any one of the calculated critical exponents. However, the extensive scaling 
of hdelity for L<^srn suggests that there may exist a subextensive contribution to hdelity, 
which is not visible at the hrst sight. To this end we take a closer look at the scaling of 
hdelity in a linear neighborhood of the multicritical point, for systems whose linear size L is 
small in the sense that 6L^1. Specihcally, we consider a neighborhood of (0, 0) quantum- 
critical point given as follows: /i=(c — 1)5, /i=(c -|- 1)5, J=(c — 1)5, J=(c + 1)5. In this 
neighborhood, the function f{k), given in (22), assumes the form 


fik) 


1 , 1_ [(c + 1)5 - p] [(c - 1)5 -p] + {(P - 1)5V _ 

2 2{(c+l)25V+[(c+l)5-p]2}V2{(c-l)252p2 + [(c-l)5-p]2}V2’ ^ ’ 


where p=cosfci -|- cos/c 2 , which after expanding to the fourth order in 5 gives 




2c 3 

+ —5^ + 
P 


5c2 + 1 


pz 


— c 



(45) 


Summing the expansion (45) over the Brillouin zone corresponding to periodic boundary 
conditions in one direction and antiperiodic ones in the orthogonal direction, with L even, 
and using the identities: l/p=0 and we obtain the following expansion 

of hdelity with respect to 6L: for 5L<Cl 


-41n%i)((0,0),5) ^a(5L)2 + 5(5L)^ (46) 

where 

a = 2-(4c^ + l)5^ b = 5c^ + l. (47) 

The expansion (46) is illustrated in Fig. 6, where the values of the coefficients a and b 
are determined numerically; it is well visible that both the coefficients are approximately 
constant in systems whose linear size satishes the inequality 5L<^1. Thus, there is no 
subextensive correction to the extensive term; the lowest order correction to the hrst term, 
proportional to is superextensive, proportional to L'^. 

Let us note here that the critical indices characterizing multicritical points of some ID 
models satisfy the condition Dv<2 but analogous analysis of the behavior of hdelity shows 
that the laws of quantum-critical scaling of hdelity are violated and one is not able to 
determine correct values of critical index v (see [16], [17]). Summarizing, it appears that 
multicritical points are not amenable to the kind of scaling analysis that we try to perform. 


5 The case of antisymmetric model 

Let us note that by applying a mean-held approximation to 

'y ^ ^l«^k.a^k,a T ^ 1 ^l^l+ej i (48) 

k,(j l,i 
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Figure 6: (Color online) The symmetric model, critical point /i=J=0, 6=fi=J. Numerical 
determination of the coefficients a and b in the expansion (46). In each panel there are 6 
plots, for three values of c: c=0 - red line, c=l - green line, c=2 - blue line, and for each c 
for three values of 6 , from left to right: (5=10“^-dotted line, 5=10“^-dashed line, 5=10“^- 
continuous line. Left panel; plots of—4 In T(i^i)((0, 0), (5)/((5L)^ a + 6((5L)^ versus L. Right 
panel: plots of (—41nT(i i)((0, 0), (5)/((5L)^ — a) /{SLY b versus L. All the plots are in 
doubly logarithmic scale. 


where Si stands for the spin operator of a spin 1/2 fermion at site I of the underlying lattice 
(for details of the notation see section 2) one obtains the Hamiltonian of the antisymmetric 
model but together with an equation relating its parameters. In distinction to our set up, 
the parameters Aj of the mean-field Hamiltonian are no longer free parameters but are given 
implicitly by those solutions of the equations 


Aj — - ai,iai+ei,-r) y 


(49) 


that minimize the ground-state energy—physical solutions in the context of system given by 
(48), where the brackets denote the Gibbs average calculated with the mean-held Hamilto¬ 
nian. When the underlying lattice is a square lattice, it turns out that the physical solutions 
satisfy the condition Ai= — A 2 , which corresponds to the so called dx 2 _y 2 pairing in theory 
of d-wave superconductivity. 

In distinction to the previously considered case of the symmetric two-dimensional model, 
where the critical points are located at straight, intersecting lines, the quantum-critical 
points of the antisymmetric two-dimensional model hll up the stripe that extends between 
the two lines |/i|=2, see the phase diagram in Fig. 7. Therefore, there are only two doubly- 
asymptotic regions of interest, where analytic formulae for G{r) can be derived, namely 
those where (/i, J)-points, with |/i|>2 and \J\ not too close to zero, approach along a p-path 
a point belonging to one of the lines \fi\=2. In these regions the asymptotic behavior of G{r) 
in the diagonal direction is 


G(r', r') 


-sgn(/i) 


J exp(-rV^(-)) 

47r^(-) r' 


(50) 


where 


&-> 


M-2 

1 + J2’ 


(51) 


that is, in the considered doubly-asymptotic region, the correlation length in the diagonal 
direction amounts to [26]- Then, in offdiagonal directions we obtained the 
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Figure 7: Phase diagram of the antisymmetric two-dimensional system in the (/i, J)-plane. 
The set of critical points constitutes the stripe |/i|<2. 


asymptotic formula 


G(ri,r2) 


Cr 


/ \ / 1 + n^ 

\ 1 -|- \ v? 


1/2 exp 


■^/^offdiag) 


cos(r0, 


(-) 

offdiag 


+ ' 




(52) 


with 


and 


1 I n" y-'V 1 1 1 \ 


q{-) ^ 

offdiag +^2 


1/2 





5 


(53) 


(54) 


provided that the points (ri,r 2 ) become remote from the origin along a ray r\jr 2 =n=const 
[26]. Sufficiently close to the lines |/i|=2, the values of \ ^2 \ \ ^2 ^ which 

determine the offdiagonal correlation length ■Cofidiag ^ direction n are given as follows (see 


[26]): 




(-) 


^ITT^+i 

1 + J 2 


Vlhl-2, 


(-) 


1 + .P 


TT 


yw-2, y-> = (55) 


1 




(-—-— 

V2 VTTj^ 

fl 1 


T'’ 

e[-\ 


(56) 

(57) 


For the considered in this section critical points 6 =jj,. In distinction to the symmetric 
model, the hierarchy of direction-dependent effective correlation lengths is not uniform in 
J [26]. However, for J>Jo, with Jo~l/4, it is the same as in the symmetric case and is 
given in (35). Consequently, for J>Jo the lower and upper effective correlation lengths and 
the associated lower and upper deviations from the critical point are given by (37) and (38), 
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Figure 8: (Color online) The antisymmetric model, critical point /i=2 and J=l, 6=fi. Plots 
of — lnT(i,o)((2,1), (5) for three values of c: c=0 - red line, c=l - green line, c=2 - blue line. 
Left panel: plots of — lnT(i_o)((2,1), 10“®) versus L. Right panel: plots of — lnT(i^o)((2,1), S) 
versus 6 for L = 10^. Both plots are in doubly logarithmic scale and the black dashed-dotted 
straight lines indicate the power-law scaling. For more details see the text. 


respectively. The critical exponents are direction independent, i^ofrdiag=i^diag=l/2, and satisfy 
the condition Du<2. In Fig. 8 we show plots of fidelity versus L and 5 in the particular case 
of J=l. Since in the scale of Fig. 8, .^iower(<5) and ^upper(<5) are quite close to each other, hence 
(Slower (-f") and Supper (-f") are close as well, one can hardly distinguish two crossover regions, so 
the mezoscopic-system regime between them is very narrow. We can say that there is just 
one crossover region whose position is given by, say, ^iower((5) (left panel) or by (5iower(L) (right 
panel). In the small-system regime of L<^iower((5) fidelity scales as which via (19) gives 
v = 1/2 in agreement with the exact value, while in the small-system regime of 5<5\o^er{L) - 
the standard 5^ scaling is observed. Then, in the macroscopic-system regime of L>^iower(5) 
the standard scaling holds, while in the macroscopic-system regime of (5>(5iower(L) fidelity 
scales as 5, which via (20) reproduces again the exact value i^=l/2. Thus, in the case of 
antisymmetric model the predictions of critical-scaling theory of fidelity match very well the 
exact results. 


6 Summary 

Hitherto, the quantum-critical-scaling theory of quantum fidelity has been verihed in a few 
ID models, which are equivalent to ID lattice-fermion gases with gauge-symmetry breaking 
interaction: an Ising chain in a transverse magnetic field [15], an anisotropic XY chain in 
a transverse magnetic held [16], and the ID version of the symmetric model considered in 
this paper [17]. It has been demonstrated that this theory holds very well in neighborhoods 
of ordinary critical points, that is those critical points that are characterized by the unique 
correlation length and the associated universal critical index. In all those cases, the quantum- 
critical-scaling theory of quantum hdelity enables one to determine the values of correlation 
lengths and critical indices. Note however that in all the cited models the condition Du<2 
is satisfied for any critical point. Only in neighborhoods of multicritical points, where the 
correlation lengths and their critical indices are not unique, breaking of the laws of quantum- 
critical-scaling theory has been observed [16],[17]. 

Our aim in this paper has been to accomplish a similar task of verihcation but for a 
higher-dimensional model. The novelty of such models, as compared to ID ones, is that 
the correlation lengths and the corresponding critical exponents u may depend on spatial 
directions and the values of u may be so large that the crucial for quantum-critical-scaling 
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theory condition Du<2 is violated-the opposite strict inequality is satisfied. In case of the 
considered here 2D quasifree lattice-fermion models, we show that this is indeed the case. To 
this end we provide analytic formulae for the large-distance asymptotic behavior of two-point 
correlation functions, the correlation lengths, and the values of the corresponding exponents 
f (a more comprehensive analysis can be found in [26]). In all those cases where we failed 
to obtain analytic results, suitable numerical results are provided. 

Some natural questions can be raised: how a multitude of correlation lengths and a 
multiplicity of critical indices at a given critical point, and or the violation of the condition 
Di/<2 is reflected in the behavior of quantum fidelity?, is it possible, in those not encountered 
in ID situations, to read off values of correlation lengths and critical indices from suitable 
plots of fidelity? 

Below, we present in a compact form the results of our attempts to answer the raised 
questions. All the data generated in our numerical calculations of fidelity lnTe(Ac,(5) are 
contained in two kinds of plots, lnTe(Ac, 5) versus system’s linear size L and lnTe(Ac, <5) 
versus a distance 5 to a critical point, in a vicinity of that point. In the paper we present 
only representative excerpts of such plots. 

The critical points of the considered models, and corresponding plots of fidelity, one can 
split into two groups: the ordinary critical points and multicritical ones. Taking for granted 
the ubiquitous dependence of quantities of interest on spatial directions, let us call a critical 
point ordinary, if for a fixed spatial direction it is characterized by the unique correlation 
length. Then, in vicinities of ordinary critical points, the examples studied in previous 
sections show that if condition Du<2 holds true in any direction, then fidelity exhibits two 
crossover regions (see Fig. 2, Fig. 4, and Fig. 8). That located at a smaller L (or 5) separates 
the small-system (quasi-critical) regime from a mezoscopic one. Its position is given by the 
lower effective correlation length .^iower((5) (or the lower effective deviation Slower (T)), which 
is the minimum over all the spatial directions of effective correlation lengths (or effective 
deviations). In the small-system regime, that is well below .Ciower(<5) (or 5iower(-h)), small- 
system scaling law (19) is satisfied. The second crossover region separates the mezoscopic 
regime from the macroscopic-system (off-critical) regime and is located by the upper effective 
correlation length Supper (<5) (or the upper effective deviation Supper (T)), which is the maximum 
over all the spatial directions of effective correlation lengths (or effective deviations). In the 
macroscopic-system regime, that is well above Supper (<5) (or Supper (-h)), macroscopic-system 
scaling law (20) holds true. If condition Du<2 is satisfied only for i/ corresponding to 
6ower(5) (or to ^upper(<5)), then Only small-system scaling (19) (macroscopic-system scaling 
(20)) is obeyed by lnTe(Ac, 5) for sufficiently large L or S but well below .^iower(<5) or (5iower(T), 
respectively (for L or S well above Supper (<5) or Supper (T)). 

After that, if in a vicinity of an ordinary critical point condition Dv<2 is violated in 
any spatial direction (see Fig. 3), then lnTe(Ac,(5) may scale, possibly in the generalized 
sense, according to a power law, and exhibit some crossovers, but neither the positions of 
crossovers are given by some ^n(<5) nor the scaling law, (19)—for sufficiently small L or 6 and 
(20)—for sufficiently large L or 6, is satisfied. In this case, of note is also the superextensive 
behaviour of — lnTe(Ac,5), which seems to defy the scaling arguments, according to which 
that should not be the case [18, 19, 20]. This curious phenomenon perhaps merits some 
further investigations into the properties of fidelity itself. 

Finally, in the case of the multicritical point of the symmetric model we have faced 
a novel situation, never before encountered in such studies: the obtained critical indices u 
satisfy the strict inequality Du>2. The behavior of In Te(Ac, S) does allow to read off neither 
correlation lengths nor critical indices (see Fig. 5). 
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